Nonequilibrium electrons in tunnel structures under high-voltage injection 
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We investigate electronic distributions in nonequilibrium tunnel junctions subject to a high voltage 
bias V under competing electron-electron and electron-phonon relaxation processes. We derive 
conditions for reaching quasi-equilibrium and show that, though the distribution can still be thermal 
for low energies where the rate of the electron-electron relaxation exceeds significantly the electron- 
phonon relaxation rate, it develops a power-law tail at energies of order of eV. In a general case 
of comparable electron-electron and electron-phonon relaxation rates, this tail leads to emission of 
high-energy phonons which carry away most of the energy pumped in by the injected current. 

PACS numbers: 73.23.-b, 74.78.-w, 74.45. +c 



Two important classes of out-of-cquilibrium phenom- 
ena in mesoscopic systems can be identified: (i) those 
that are described by quasi-equilibrium quasiparticle dis- 
tributions where electrons and phonons have well de- 
fined, although in general different, temperatures and 
(ii) those where quasiparticle distributions deviate sig- 
nificantly from their equilibrium form and where the no- 
tion of temperature cannot be introduced. The correct 
interpretation of the data observed in the particular ex- 
periment requires the proper attribution of the specific 
phenomenon to either of these classes. Superconduct- 
ing mesoscopic structures such as tunnel and Joseph- 
son junctions arc among those that are being investi- 
gated most intensively. The common tools for infer- 
ring information about nonequilibrium states in such 
systems are studies of the electron-phonon interaction 
rates [H, @, 3 , IE), |j| and of the electron-phonon energy 



relaxation 




1121 . 113l | . A standard setup is 
an island connected to leads (they all can be supercon- 
ducting or normal) via tunnel contacts and driven out 
of equilibrium by strong electron injection under bias 
voltages V such that eV significantly exceeds both the 
temperature and the superconducting gap A. In many 



experiments (see, e.g., Ref. Il3l). a quasi-equilibrium dis 



tribution is created with an electronic temperature T c 
essentially higher than the bath temperature T p h kept 
low by efficient cooling. These temperatures determine 
the energy transfer from electrons to the phonon bath. 
The quasi-equilibrium establishes when a big difference 



between fast electron-electron, 



7e 



and slow electron- 



phonon, 7 c _ph, relaxation rates exists at energies e ~ T c . 
These conditions arc usually fulfilled in aluminum sam- 
ples at sub-Kelvin temperatures. 

The more common is the situation, however, where, 
in contrast to Al samples, the ratio 7 c -o/7o-ph at such 
temperatures is not very high (as illustrated, e.g., by 
the data in Ref. In this Letter we investigate for- 



mation of the electronic distribution in the general case 
and show that for moderate ratios 7 e -o/7c-ph the condi- 
tions of quasi-equilibrium can be easily violated. The 
distribution formed under a high-voltage injection can 
be characterized by an effective electronic temperature 
T c > Tph only within a low-energy region around e ~ T c 
provided 7 e - e (T e ) is much larger than 7 c - P h(T ). We find 
that at energies of order of eV ^> T c , the distribution 
has a long power-law tail, which crucially changes the 
transport properties of the entire electronic system hav- 
ing much lower effective temperature T e . Since the tail 
energies exceed A, our results are general and apply to 
both normal and superconducting junctions. While the 
derived behavior somewhat resembles the well known 
electron runaway in semiconductors [l4| , it is drastically 
different from such nonequilibrium effects as supercon- 
ductivity stimulation [IB}, nonequilibrium proximity ef- 
fect in Josephson junction (see, e.g., [l6| and references 
therein), etc., where the deviation from equilibrium is 
maximal right in the energy range near A. 

We consider the energy exchange between the electron 
and phonon subsystems relevant to a rich variety of ex- 
periments @, i, \& E El, E El- Since 7 e . ph (e) is well 
known to grow faster than 7 e _ e (e) with increasing energy, 
there exists certain energy, e*, at which these two rates 
match. One would expect that nonequilibrium effects in 
the energy transfer are small as long as T c <C e* . We show 
that it is indeed the case if eV <C e* and the electron- 
electron interaction dominates in the entire nonequilib- 
rium region. However, at eV 3> e* a crossover from the 
electron-electron to the electron-phonon mechanism of 
relaxation takes place as a function of energy. In this case 
the power law tail in the electron distribution leads to 
emission of high-energy phonons, which carry away most 
of the energy provided by the injected current. Accord- 
ingly, the energy emitted via thermal phonons becomes 
a much smaller fraction of the total inserted energy, im- 



2 



plying that interpretation of the experimental data based 
on the quasi-equilibrium distribution function with some 
effective temperature does not apply. 

Setup. Having in mind mostly superconducting de- 
vices, we consider a junction consisting of a supercon- 
ducting (or normal) macroscopic-size island (i) connected 
via small-area high-resistance tunnel contacts with two 
superconducting (or normal) leads L\ and Li. The bias 
voltage V is assumed high in the scale of characteristic 
energies of the particular experiment. For the system 
mentioned above, the proper inequality is eV 3> T e ,A. 
For the energies T e ,A <C e < eV, which we are inter- 
ested in, the normal-state equations for the electron- 
phonon, electron-electron interactions, and the density 
of states can be used both for a normal and super- 
conducting junction. The electron temperature is as- 
sumed uniform along the sample. We use the clean- 
limit approximations for 7 _ c (e) and 7 e - P h(e) because im- 
purity scattering does not significantly rcnormalize the 
electron-electron or electron-phonon interactions for en- 
ergies of interest. The influence of the impurity scatter- 
ing on the electron-phonon relaxation is controlled by 
the parameter q£ where £ is the electronic mean free 
path and q = e/Hs is the wave vector of an emitted 
phonon with energy e 17, [3, Taking sound ve- 

locity s ~ 5000 m/s and £ ~ 20 nm as in Al sam- 
ples of [H we get q£ ~ (0.5 K" 1 ) e/k B . Therefore, at 
e/ks ~ eV/ks » 1 K the clean limit is appropriate. 
The impurity-induced renormalization of the electron- 
electron interaction [2(| becomes important for energies 
e < H 3 /t(pf£) 2 ~ £F(ft/PF^) 3 and can also be neglected 
for e ~ eV. Under these assumptions we can derive a 
linear equation for the distribution function and solve it 
exactly for the most relevant situations. 

We consider a symmetric structure with voltages at the 
leads Vl 1 = —Vl 2 = V/2 such that the chemical potential 
of the island \ii = by symmetry, while for the leads 
Hl 1 = —hl 2 = Mi — —eV/2. It is convenient to write 
kinetic equations for odd and even components of the 
distribution function n(e) defined as /i(e) = n{— e) — n(e), 
and /2(e) = 1 — n(e) — n(—e). If the leads are in thermal 
equilibrium at temperature Tj,, we have f[ Ll%> = f[ L2 ^ = 



f[ L) and f™ 



AL) 
Jl,2 



tanh ■ 



(£2) 



f^ L ' where 
^itanh^^ 2 



2Tr, 



2TY 



In what follows we assume that the leads and the phonon 
bath are at zero temperature, Tl = T p h = 0, and that all 
the emitted phonons arc immediately removed from the 
sample due to ideal heat contact to the substrate. 

One can check that the even component of the distri- 
bution in the island vanishes by symmetry, /2(e) = . 
The kinetic equation for the odd component is 



r(ph) 



where ', j} ', and j[ ph ^ arc, respectively, the tunnel, 
electron-electron, and electron-phonon collision integrals 
in the island. The latter describe relaxation of the distri- 
bution driven out of equilibrium by the tunnel source 21 



(2) 



that contains the distribution in the leads, f[ L \ and 
in the island, f\. The (identical) tunneling contacts 
are characterized by an effective tunneling rate 77 = 
(4^e 2 f2i?) -1 where v = v{Ef) is the normal density of 
states in the island, Q is its volume, and R is the contact 
resistance. In what follows we consider high contact re- 
sistances, i. e., small r\ (the estimate will be given later). 

Electron-phonon relaxation. For small 77, the distri- 
bution function determined by Eq. (fTJ) is close to the 
thermal, /i(e) « tanh(e/2T e ), with a certain electronic 
temperature T . For energies e ~ T c the deviation from 
quasi-equilibrium is negligible. However, at e 3> T c where 
tanh(e/2T e ) ps sign(e), the injection-induced deviation 
becomes essential. We put /i(e) = tanh(e/2T ) + 5 fie) 
where |<5/| <C 1. At T p h = 0, the electron-phonon colli- 
sion integral [22| in the island vanishes for f\ (e) = sign(e). 
Therefore, neglecting exponentially small terms we get 
for energies |e|, |e + cj| 3> T c and e > (see also Ref. [23[) 



(ph) _ 7e-ph(r c ) 



e 3 Sf(e) 



duui 2 5f(e + oj) 



(3) 



where 7 e - P h(Tc) = ^K- P hT 3 /2h{spF) 2 is the electron- 
phonon relaxation rate at the electronic temperature 
T , while A e -ph is the interaction constant. We use T e 
simply as a convenient energy scale. In fact, the ratio 
7e- p h (T e )/T c 3 is independent of T e . 

Electron- electron interaction. The electron-electron 
collision integral satisfies the energy conservation 
law, ej[ e \e) de = 22 1. For large energies 

lei, |ei|, leal, 1^3 1 ^ T Q and e > it has the form 



7c-c( r c) 



2 Sf(e) 



— 3 / duo uj Sf(e 




■ u>) 



(4) 







(1) 



Here 7 C - C (2c) = ttA c _ c T c 2 /8hEp is the electron-electron 
relaxation rate at T e , A - is the interaction constant. 
The ratio 7 e - e (T , e )/T e 2 is independent of T e . 

Distribution function. The electron-electron relax- 
ation dominates at very low energies, where the distri- 
bution has a thermal form with an electronic temper- 
ature T e . At higher energies, a deviation from thermal 
behavior develops due to the reduced role of the electron- 
electron interaction. Consider e > 0. Since /i(e) w 1 
for £>r e the tunnel collision integral Eq. ((2]) becomes 
= —4?/ 9 (eV/2 — e). This form of the injection term 
suggests that at e > T e , Sf(e) = -0(e) 6 (eV/2 - e). 
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Using Eqs. d2j) — (J4j) , the kinetic equation ([T]) becomes 



1 



eV/2- 



dui uj 2 (j)(e + uj) 



eV/2-e 



du) UJ <p(e + lu) 



1.(5) 



Here e 3 ph = 4 7? T c 3 / 7c . ph (T c 



and 



Equation (J5]) has a characteristic energy scale 



4r/T 2 / 7c _ c (T e ). 



fc ph 



Ae ~ r [7 c _ (T )/7 e _p h (T e )] 



(6) 



such that 7e- P h(e*) = 7c-c(e*)- Therefore, the electron- 
phonon interaction dominates for e 3> e* while the 
electron-electron interaction dominates for e <C e* • As it 
is explained earlier, we consider here the situation when 
7e-e(Te) > 7o- P h(T' ), thus e* > T c . 

The integral equation ([5]) can be transformed into a dif- 
ferential one by triple differentiation over energy, which 
can then be easily analyzed and solved numerically The 
boundary conditions are obtained by putting e = eV/2 at 
each step. In the situations of dominant electron-electron 
relaxation, 7 C - P h — * (i.e., e p h — ► oo), it can be reduced 
to the second-order differential equation for a function (j) e 
determined entirely by the electron-electron interaction 



de 



[eVe(e)] =60 e (e): 



(7) 



with the boundary conditions e 2 <p c (e) — 2e 2 and 
{d/de)[e 2 (t) c {e)] = at e = eV/2. Solution of Eq. © 
with these boundary conditions is 



4e 2 
5 



3 (eU/2) 2 
2 e 4 



1 (eU/2) 3 
Fast electron-phonon relaxation, 7e - p h — * oo, leads to 
d 3 



(8) 



de 3 



[e 3 </) ph (e)] +60 ph (e 







(9) 



3e 3 h , (d/de)[e 3 ph ( e )] 



with the conditions e 3 p h(e) 

(cP/efe 2 )[e 3 c/> ph (e)] = at e = eV/2. The solution is 
18 



P h(e) = -(e ph /eyT(eV/2e) , 



(10) 



- cos(v / 21na;) + sin(v / 21na;) 
6 3 



It was obtained in [23j and used in |24( for studies of 
nonequilibrium created by absorption of a high-energy 
photon. 

At comparatively low voltages, eV -C e*, the electron- 
electron interaction dominates in the entire energy range 
< e < eV/2 and the distribution function obeys Eq. ((8]). 
At high voltages, eV 3> e*, one can discriminate between 
two regions with different relaxation mechanisms with a 



crossover between them at e ~ e*. For < e <C e* the 
electron-electron interaction dominates. For e* <C e < 
eV/2 the electron-phonon interaction wins, and the dis- 
tribution function approaches Eq. (fl"0|) . In both cases the 
electronic distribution has a long power-law nonequilib- 
rium tail 4>{e) ~ (e/e) 4 at energies T e -C e <C eV. For 
eU < e*, one has e = (eeeU) 1 / 2 , while e = (ej^eV) 1 / 4 for 
eV 3> e*. The deviation from equilibrium becomes of the 
order unity for e < e. Thus the low energy distribution 
can be thermal only if e -C T c . This condition reads 



eU/e* « £ • max{£, 1}, $ = T 2 /e c e* 



(11) 



Therefore, one can interpret the data inferred from the 
experiment in terms of a quasi-equilibrium electronic 
temperature T e only if the condition (fTTjl is fulfilled. Oth- 
erwise, the distribution is not thermal even at small e; 
instead, the scale of its variation, and thus the apparent 
"temperature" , is determined by e. 

Energy balance. Even if the quasi-equilibrium condi- 
tion (fTT|) holds, the long power-law tail in the distri- 
bution can strongly influence electronic processes at a 
much lower temperature T c . An important example is 
the electron-phonon energy relaxation. Using Eq. ([I]) one 
finds the energy balance per unit volume of the island 



2v 



i-eV/2 reV/2 

/ e Jj de + 2v / eJ^ ph) de = 
Jo Jo 



(12) 



The electron-electron collision integral vanishes due to 
the energy conservation. With Eq. Q this gives 



P eq (T e ) 



Pnc q = PV 



(13) 



Equation (|13[) determines T e as a function of the in- 
jected power Py = U 2 /4i?f2. The latter is obtained from 
the first integral in Eq. (|12p under the assumption that 
eV T c . This is the power (per unit volume) deposited 
into the island in a setup with two contacts of total resis- 
tance 2R under the total voltage bias V. It is half of the 
total input energy, the other half of which goes into the 
leads. The l.h.s. of Eq. (fT3|) comes from the second inte- 
gral in Eq. (fT2|) and consists of the energy transfers, P cq 
and -Pneq, to zero-temperature phonon bath by thermal- 
ized and by nonequilibrium electrons, respectively. The 
energy transferred by thermal electrons, P cq (T ), is deter- 
mined by the electron-phonon interaction integrated over 
the low-energy domain e > T e with f\ = tanh(e/2T e ) ne- 
glecting the nonequilibrium correction. For a clean nor- 
mal metal, P oq (T ) = ET e 5 where S w 78is\/{sp F ) 2 h Jl. 
For clean superconductors, P eq (T e ) was calculated in [1 31 ] - 
The energy transferred by non-thermal electrons is 



P 



ncq 



T3 - 



j-eV/2 




/ ede 




'o 


3 Jo 



eV/2-e 



dui UJ 4>(e + lu) 



(14) 



The main contribution here comes from e ~ eV. 
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For voltages eV <C e* when the electron-electron inter- 
action always dominates, the non-thermal contribution 
can be calculated using Eq. ([5]). The e -4 term in Eq. 
(dU) is cancelled out leading to P ncq = (eV/3e*)Py, and 
Eq. (0 yields 

P eq (T e ) = P v [l-eV/3e*} . (15) 

The second term in the brackets is a small correction: Al- 
most all the injected energy is absorbed and then trans- 
mitted to the phonon bath by thermal electrons. 

The situation is totally different at higher bias volt- 
ages, eV 3> e*. In this case almost all the injected power 
is absorbed by high energy electrons, and T e is less sensi- 
tive to V. Using p h(e) from Eq. (fTO|) . which satisfies the 
kinetic equation without the electron-electron collision 
integral, the energy balance Eq. (fl"3|) becomes 

P cq (T c ) = ^ / e de [0 ph (e)-0(e)] 

£ ph JO \ J 

reV/2-e \ 

-J uj 2 duj [0 p h(e + w) - (f>(t + w)] j . 

The function 0(e) satisfies the full kinetic equation {5]). 
The main contribution to the integral comes from e ~ 
eV 3> e* where Eq. (0 coincides with that for small 
electron-electron interaction within the accuracy e* /eV. 
Therefore (0 p h—0)/0ph ~ e*/eU. This results in a slower 
dependence of P eq on the bias voltage, 

P e < l (T e )~P v (e*/eV)=vr l e*eV. (16) 

Expressions (fT5"|) and (fT6")) are the main result of our Let- 
ter. They reveal a crossover from oc V 2 to oc V behavior 
that occurs at eV ~ e*. This is confirmed by numerical 
solution of the kinetic equation (|5|) . 

Discussion. The data for aluminum samples [3, [Til . 
El suggest 7o (T c ) - 10 8 it 1 , 7 P h(T c ) ~ 10 6 s"\ while 
?; = 10 s _1 for the samples of Ref. [l3j|. Therefore, 
e*/T c ~ 10 2 . For voltages eV ~ 10 2 T c used in 0, 
the electron-electron interaction always dominates, while 
the deviation from equilibrium remains small. Thus the 
quasi-cquilibrium conditions are well satisfied. 

However, the above situation is rather an exception 
to the rule - it is in a sense unique to aluminum sam- 
ples. In other materials the electron-phonon relaxation 
rates are normally much larger [l[. Consequently, their 
ratios to the corresponding electron-electron relaxation 
rates are not as low as in Al. Therefore, the crossover 
from electron-electron to electron-phonon dominated re- 
laxation occurs at much lower energies, and should be 
observable in experiments. Although the deviation from 
the thermal distribution could be still small within the 
low-energy range e ~ T e , the high-energy tail of the dis- 
tribution at e ~ eV becomes essential for the energy 
exchange between the electron and phonon subsystems. 



The injected power will then be mostly absorbed by high- 
energy phonons with energy T c <C e < eV. 

To summarize, we have established the quasi- 
cquilibrium condition Eq. pip , and demonstrated that 
the energy exchange rate P eq (T e ) between thermal elec- 
trons and the phonon bath experiences a crossover from 
Eq. (fTS")) to Eq. (fTrJ)) as a function of the bias voltage. 
The crossover energy determines the relative strength of 
the electron-electron and electron-phonon interactions. 
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